Application of the generalized two-center cluster model to 10 Be 
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A generalized two-center cluster model (GTCM), including various partitions of the valence nu- 
cleons around two a-cores, is proposed for studies on the exotic cluster structures of Be isotopes. 
This model is applied to the 10 Be=a + a + n + n system and the adiabatic energy surfaces for 
a-a distances are calculated. It is found that this model naturally describes the formation of the 
molecular orbitals as well as that of asymptotic cluster states depending on their relative distance. 
In the negative-parity state, a new type of the a-(- 6 He cluster structure is also predicted. 
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PACS numbers: 21.60.Gx,24.10.Eq,25.60. Je 

Recent experiments by Freer et ai., revealed the 
existence of the interesting resonant states, which dom- 
inantly decay to 6 He s . s .+ 6 He g . s . and a+ 8 He g . s . chan- 
nels in the 10 to 25 MeV excitation energy interval of 
12 Be. Similar resonant states, decaying to He-isotopes 
such as 6 He and 8 He, have also been observed in other 
Be isotopes of 10 Be Q and 14 Be 0. These experiments 
strongly suggest the existence of exotic cluster states con- 
sisting of the respective He-isotope clusters. 

Theoretically, Kanada-En'yo and her collaborators Q 
have studied a wide range of light nuclei within the anti- 
symmetrized molecular dynamics (AMD), and predicted 
many kinds of cluster structures in isotope chains of the 
elements Be, B and C. In particular, much attention has 
been concentrated on Be isotopes in which a motion of 
valence neutrons couples to the two-a structure of 8 Be. 
Based on the molecular orbital model (MOM), Itagaki 
et al. has extensively studied the low-lying states of Be 
isotopes to clarify the roles of valence neutrons on a de- 
velopment of clusterization. The molecular orbitals, such 
as 7T and a orbitals associated with the covalent binding 
of atomic molecules, have been shown to give a good 
description for the low-lying states of Be-isotopes as em- 
phasized in the AMD studies 0,0. 

On the other hand, Descouvemont and Baye 
have applied a traditional cluster model ( 6 He+ 6 He) + 
(a+ 8 He), assuming the substructures of 6 He s . s , and 
8 He s . s ., to illustrate the resonances observed in 12 Be. 
However, it is plausible to assume 6 He and 8 He nuclei to 
be stable and frozen clusters because the energy for two- 
neutrons separation is quite small in both nuclei. Thus, 
it is properly considered that the valence neutron of 6 He 
and 8 He are associated with companion clusters when 
they approach each other, which is similar to the cova- 
lent binding of valence neutrons in the low-lying states. 

The purpose of this letter is to study the molecular 
structures in low-lying and high-lying states in a unified 
way by taking into account couplings between a valence 
neutron's motion and a relative motion of the a clus- 
ters. To achieve this purpose, we propose a new frame- 
work of a generalized two-center cluster model (GTCM) 
where we can describe atomic-orbital motions of valence 



neutrons around individual a clusters on the same foot- 
ing with molecular orbitals, being a single-particle mo- 
tion around two a-cores. We apply this model to the 
w Hc=a+a+n+n and discuss its applicability to stud- 
ies of the molecular structure in both the low-lying and 
high-lying states. 

The basis functions of GTCM for 10 Be are given as 

( 10 Be; S) 
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where C^^(S) means the product of d^ K • d^ K - The 
a-cluster wave function ipi(a) (i=L,R) is given by the 
(0s) 4 configuration in the harmonic oscillator (HO) po- 
tential with the relative distance-parameter S. The po- 
sition of an a-cluster is explicitly specified as the left 
(L) or right (R) side. A single-particle state for valence 
neutrons around one of a clusters is given by an atomic 
orbitals, ip(i,p n , r) with the subscripts of a center i (=L 
or R), a direction p n [n=x, y, z) of Op-orbitals and a 
neutron spin r {—] or {). In Eq. the index m(n) is 
an abbreviation of the atomic orbital {i,Pm t). The basis 
function $^ n ( 10 Be; S) with the parity 7r is projected to 
the eigenstate of the total spin J and its intrinsic angu- 
lar projection K by the projection operator Pj^. Various 
linear combinations of ^(S) can be shown to repro- 
duce not only molecular orbital configurations but also 
cluster-model states of 4 He+ 6 He and 5 He+ 5 He. 

First, we illustrate that the molecular-orbital configu- 
ration of (c + ) 2 can be constructed from a linear combi- 
nation of the basis function for instance. Since a + orbital 
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is expressed as ip(L,p z ,r) — ip(R,p z ,r), the wave func- 
tion of 10 Be with (cr + ) 2 of valence neutrons are written 
as follows: 

$+( 10 Be; S) = A{il> L (<x)M<x) 

x(v(L,Pz, t) - <p(R,Pz, T)) • (v(L,p z , I) - <p(R,Pz, I))} 

= A{ip L (a)(p(L,p z ,1)(p(L,p z ,i) ■ ip R (a) 
-ip L (a)ip(L,p z ,1) ■ i/j R (a)(p(R,p z ,l) 
-ip L (a)ip(L,p z ,l) ■ i/j R (a)ip(R,p z ,1;) 
+Ma) ■ i/>R(a)tp(R,p„1)<p(R,p t ,l)}- (3) 

This expression means that a 10 Be wave function with the 
molecular orbitals can be described by a linear combina- 
tion, such as ( 6 He-a)+( 5 He- 5 He)+(a- 6 He), because the 
clusters of 5 He(=a+n) and 6 He(=a-t-2rt) are given by 
A{tpi(a)ip(i,p n ,T)} and A{i(Ji(a)<p(i,p n ,T)tp(i,p n , ,t')}, 
respectively. This linear combination can be exactly con- 
structed from basis functions of <3E>^ „( 10 Be; S). 

Next, we show that the basis functions of Eq. @ 
can also describe the cluster-model states in which 
neutron's orbitals have a definite spin around one 
of a-cores. For instance, 0p 3 / 2 jz=+3/2 is ex- 
pressed as {ip(i,p x ,1) — iip(i,p y ,1)}. Thus, the 3/ 
2~ states of 5 He with the 0p 3 / 2 neutron is writ- 
ten as Y^n-A-iipiiot) ■d n cp(n)}= A{ipi(a) ■ ip(i,p x ,1)} - 
A{tpi(a) ■ i(f(i,p y , T)}. Similarly, the 6 He(/ 7r ) clus- 
ters with a definite intrinsic spin-parity I n , such as 
[(0p3/2) 2 ]o+ j can be constructed from a certain linear 
combination of J2 m n A{4>i(a) ■ d m ip(m) d n ip(n)}. There- 
fore, the wave function of Eq. can describe the cluster- 
model states of [a ® 6 Hc(/ 7r )] and ^He^ 1 ) ® 5 He(Jj 2 )]. 

The wave function of 10 Be is finally given by taking a 
superposition over the relative distance-parameter S and 
the intrinsic angular projection K as 



^'( 10 Be) = / dS X! *f ( 10 Be;5) 

= fdSj2cf K (S)<S> J K(S) (4) 

J K/3 

with /? = (m, n). The coefficients Cp K (S) are deter- 
mined by solving a coupled channel GCM (Generator 
Coordinate Method) equation 0: 

f dsJ2cf K (s) 

x ($£ K '(S') \H - E J * | $£K(S))=Q . (5) 

To see the coupling properties, we solve Eq. JSJl in a 
step by step. First, we solve Eq. JSJ at a fixed S. Namely, 
we solve 



0K 

x (*£ K \S) | H - E J '(S) | <S>f K (S))=0 .(6) 

The eigenvalue E J (S) is a function of the relative 
distance-parameter S, and then we call the solutions of 
energies and wave functions "adiabatic energy surfaces" 
and "adiabatic eigenstates" , respectively. The calculated 
adiabatic energy surfaces for the J 7r =0 + state are shown 
by open circles in Fig. ^ Here, for the nucleon-nucleon 
interaction, we adopted the Volkov No. 2 with the Majo- 
rana parameter m=0.576 and without the Bartlett and 
Heisenberg exchanges. Due to the reduction of the ma- 
jorana parameter, the total binding energy is gained and 
it becomes easy to see the continuum states. We also 
employed the G3RS interactions for the spin-orbit parts. 
The radius parameter b of HO wave functions for a clus- 
ters and valence neutrons is taken as 1.44 fm. 

At the asymptotic distance (S— >oo), where two 
a-cores are completely separated, we can define 
the asymptotic channels such as [ 4 He+ 6 He(/)]i and 
[ 5 He(/i)+ 5 Hc(/2)]/L, in which individual clusters have 
intrinsic spins (Ii, I2) and coupled with the channel spin 
/ (I=Ii + 12) and the relative one L. We call the cou- 
pling scheme of these asymptotic channels as "a cluster- 
coupling scheme" . The solid and dotted curves shown in 
the right part of Fig. are the expectation values < H > 
of the [ 4 He-(- 6 He(/)] L and [ 5 He(/i)+ 5 He(/ 2 )] /L cluster- 
coupling schemes, respectively 

From Fig. ^ we can see that in S > 6 fm, the calcu- 
lated energy surfaces are completely the same as those 
of cluster-coupling wave functions. This means that, in 
this region, the valence neutrons are localized at one of 
a-cores and rotation of two clusters is de-coupled to each 
other. At the asymptotic region, therefore, each nucleus 
keeps its isolated states and weakly coupled to each other. 
On the other hand, in S < 6 fm, the energies for the 
cluster-coupling scheme deviate from the adiabatic en- 
ergy surfaces. 

We study the adiabatic eigenstates in an internal re- 
gion (S < 4 fm) in the view of the molecular orbital 
formation. As shown in an example of (c + ) 2 in Eq. J3J), 
the wave function of molecular orbitals is expressed by a 
linear combination of different kinds of cluster- wave func- 
tions, where each cluster has no good angular momenta. 
In the molecular orbitals, the valence neutrons are mov- 
ing around two a-cores with a specific direction in respect 
to the a— a axis. Such a configuration of the system is 
called as "a strong-coupling scheme" . The overlap of the 
adiabatic eigenstates with the wave functions of molec- 
ular orbitals identifies the dominant components in the 
adiabatic eigenstates as shown in Fig. The adiabatic 
eigenstates connected by the thin-solid curves at the in- 
ternal region have the common dominant-components of 
various kinds of molecular orbitals. 
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FIG. 1: Dynamical transition from the molecular orbitals to 
the cluster-coupling states in the adiabatic energy surfaces. In 
the right part, the solid and dotted curves show the energies 
of the a+ 6 He and 5 He+ 5 He cluster states, respectively. In 
the left part, the surfaces of A, B, C and D have a dominant 
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FIG. 2: The same as Fig. but for the negative par- 
ity states (J v =l~). In the internal region of S <4.5 
fm, the lowest, second, third and fourth curves connected 
by a thin-solid line show the surfaces having a dominant 
component of (>3~ /2 o"+ /2 )if = i, Or /2 °"^ /2 )if=i, (tt~ /2 <t+ 2 ) k=0 , 
( 7r 3/2 7r i/2) K=0 ' respectively. The surface with a thick curve 
has a dominant component of the [a+ 6 He(Oj")]L=i channel. 



TABLE I: Dominant component in the eigenstates along to 
the lowest adiabatic surface A of Fig. In up-most row, 
the distances of the calculated eigenstates are shown. The 
dominant basis-state is shown in the second row, while its 
population is done in the lowest one. The distance parameter 
of 5=6 fm corresponds to the "intermediate-coupling region" . 
See text for details. 



S (fm) 


2 


4 


6 8 


Main component 
Squired ampli. 


0.93 


( 7T 3/ 2 y 

0.81 


Mixed [a+ b He(0+)]i,=o 
0.99 



Table [I] shows the dominant molecular-orbital config- 
urations in the adiabatic eigenstates along to the lowest 
surface A of Fig. ^ The squared amplitudes listed in Ta- 
ble [I] have the meaning of the probability of finding the 
system in a molecular orbital, because an orthogonality 
among the different molecular orbitals is good in spite of 
the anti-symmetrization effect. This is due to the fact 
that the molecular orbitals are similar to the deformed 
shell-model orbitals in the internal region. 

We can see that the eigenstate has a dominant compo- 
nent of (7i"^/ 2 ) 2 a ^ S=2 and 4 fm. On the other hand, 
at the asymptotic region (5 ~ 8 fm), the molecular 
orbitals are strongly mixed with each other, and the 
[a+ 6 He(0^)]L = o channel dominates. Thus, the cluster- 
coupling scheme becomes a good basis state for de- 
scribing the system at the external region, while the 
strong-coupling one at the internal region. At an in- 
termediate region (S'=5~6 fm), we find that the dom- 
inant component of the lowest surface A is changed from 



( 7r 3/2) 2 to (c^) 2 At this distance, an am- 

plitude of the [cv+ 6 Hc(0^)]l=o channel in the cluster- 
coupling schemes amounts to 61%. This means that, 
in .5=5^6 fm, the eigenstate has "an intermediate cou- 
pling scheme" between the strong-coupling scheme and 
the cluster-coupling one. 

Finally, the GCM equation (JSJ is solved by employing 
the basis states ranging from S=l fm to S=9 fm with the 
mesh of 0.5 fm. We obtain the result that the lowest three 
GCM solutions have energy gains of about 1^2 MeV, 
and the dominant amplitudes around the respective local 
minimums in the adiabatic energy surfaces A, B and C. 
Therefore, these solutions of the 0^~, Oj" and 0% states are 
concluded to have the molecular-orbitals configuration of 
( 7r 3/2) 2 ' ( CT i/2) 2 an< ^ ^ 7r i/2) 2 ' res P ec tively Q. 

Similarly, we calculate the J n ~l^ state. The calcu- 
lated energy surfaces are shown in Fig. We focus on 
the lowest five surfaces and investigated their intrinsic 
structure, because there is no definite local-minimums in 
higher surfaces. As shown by two vertical dotted-lines in 
Fig. [21 we find the dynamical transition from the strong- 
coupling scheme to the cluster-coupling one in the adi- 
abatic surfaces. However, the first excited surface con- 
nected by a thick curve has an almost pure-component of 
the [a+ 6 He(0^)]L = i channel. Thus, this excited surface 
is made by the pure cluster-coupling state of this channel 
in a wide range of the a-a distance. 

The energy gain of the lowest three solution due to 
GCM is about 1~3 MeV depending on the states. In the 
lowest and the third l - states, we find that the wave 
functions distribute around the respective local mini- 
mums of the adiabatic surfaces. Thus, the intrinsic struc- 
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FIG. 3: Band structures in Be. The bands shown by 
the double circles and the white squares have (7r.T, 2 ) 2 and 
( 7r 3/2 (J i/2)' respectively, while that with the white circles does 
(crjy 2 ) 2 . The solid symbols have a developed a+ 6 He cluster 
structure. The inset is a same figure but for the observed 
states 0,0- The spins for the states in parentheses are ten- 
tatively suggested in Ref. Q. 

tures of the lj~ and states are explained with the 
molecular orbitals of {^p<Jip)K=i and (Tr~^ 2 a^ 2 )K=i, 
respectively. The 1 2 state has the main components of 
the first excited surface in Fig. which is interpreted in 
terms of the cluster-coupling states of [a+ 6 He(0f ))l=i- 
To see band structures in 10 Be, we solve GCM Eq. (|5J) 
for the higher spins with a natural parity (— 1) . The 
calculated bands are shown in Fig. |3] with an inset show- 
ing the respective observed states [ij, Q • The moment of 
inertia of individual bands is well reproduced. Further- 
more, our model predicts the existence of the higher spin 
states which were suggested by a recent experiment 
In the (cjy 2 ) 2 band shown by the white circles, we find 

the enhancement of the [a+ 6 He(2^)]£=4 state at a maxi- 
mum spin (solid circle) . In the present calculation, all the 
excitation energies is higher than those of the observed 
states. This can be improved by optimizing the nucleon- 
nucleon force and the radius b to reproduce the cluster's 
threshold energies. The reproduction of threshold ener- 
gies is quite important for the study on the decay width. 
The quantitative analysis including the decay width will 
be given in forthcoming papers. 



In summary, we proposed a generalized two-center 
cluster model (GTCM) and discussed its application to 
the 10 Be—a+a+n+n system with J n =0 + and 1~. The 
adiabatic energy surfaces depending on the a— a dis- 
tance were calculated. The adiabatic eigenstates have 
the molecular orbital configuration at an internal region, 
while they becomes the cluster-coupling states at an ex- 
ternal one. The middle distances correspond to the tran- 
sitional region having an intermediate coupling-scheme 
between the internal regions and the external ones. It 
should be noticed that, in our model, both schemes of 
the molecular-orbitals and cluster-coupling schemes can 
be naturally described in an equal footing without any 
difficulty relevant to the double-projection procedure. 

Finally, we solved the coupling between the relative 
motions of clusters and the intrinsic motion of valence 
neutrons. The low-lying + states and the lowest 1~ 
one have the respective molecular-orbital configurations, 
which are consistent with previous studies based on 
MOM 4] and AMD 0. In addition, our model pre- 
dicts the possible appearance of the [a+ 6 Hc(2 ] f )}l=a 
cluster-structures in the lowest 6 + state. Furthermore, 
we theoretically obtain the second 1~ state with the 
[a+ 6 He(0f )]l=i structure and expect to observe it ex- 
perimentally. To see correspondence with experiments, 
it is necessary to analyze resonant states above threshold 
energies. We are now going to tackle this problem. 

In conclusion, we can say that GTCM well describes 
the low-lying molecular orbitals obtained by other the- 
oretical models 0, 13 • Furthermore, this model nat- 
urally reproduce the asymptotic cluster-coupling states 
which are described by the traditional cluster model 6] in 
an equal footing. These results indicate that the present 
GTCM is applicable to the study of low-lying and high- 
lying states in 10 Be. Since the a+a+2n model for 10 Be 
is a special case (x=2) of a general a+a+xn model, it is 
very easy to apply GTCM to systematic studies of res- 
onances observed in excited Be-isotopes as well as their 
low-lying states. In particular, its application to 12 Be 
= a+a+in is very interestingbecause of their accumu- 
lated experimental results 0,13- A direct extension to 
the a+a+4n model for 12 Be is is now under progress. 

The authors would like to thank N. Itagaki, A. Ohnishi 
and Y. Sakuragi for their valuable discussions. One of 
the present authors (M. I. ) also would like to thank the 
Japan Society for the Promotion of Science (JSPS) for 
financial support. This work was performed as a part of 
the "Research Project for Study of Unstable Nuclei from 
Nuclear Cluster Aspects " at RIKEN. 



[1] M. Freer et al, Phys. Rev. Lett. 82, 1383 (1999); Phys. 

Rev. C 63, 034301 (2001), and references therein. 
[2] A. Saito et al, Suppl. Prog. Theor. Phys., 146, 615 (2003). 
[3] See for examples, Y. Kanada-En'yo, H. Horiuchi and A. 

Ono, Phys. Rev. C 52, 628 (1995) ; Y. Kanada-En'yo and 



H. Horiuchi Phys. Rev. C 55, 2860 (1997), and references 
therein. 

[4] N. Itagaki and S. Okabe, Phys. Rev. C 61, 044306 (2000) 
; N. Itagaki, S. Okabe and K. Ikeda, Phys. Rev. C 62, 
034301 (2000), and references therein. 



■5 



[5] Y. Kanada-En'yo, H. Horiuchi and A. Dote Phys. Rev. C 
60, 064304 (1999). 

[6] P. Descouvemont and D. Baye, Phys. Lett. B505, 71 
(2001) ; D. Baye, P. Descouvemont and R. Kamouni, Few- 
Body Systems 29, 131 (2000). 



[7] H. Horiuchi et al, Suppl. Prog. Theor. Phys. 62, 1 (1977) 

and references therein. 
[8] Y. Ogawa et al, Nucl. Phys. A673, 122 (2000). 
[9] J. A. Liendo et al, Phys. Rev. C 65, 034317 (2002). 



